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ABSTRACT 

An asymptotic approach to solution of the inverse problems of 
remote sensing is presented. It consists in changing integral opera- 
tors characteristic of outgoing radiation into their asymptotic 
analogues. Such approach does not add new principal incertainties 
into the problem and significantly reduces computation time that 
allows to develop the real (or about) time algorithms for interpre- 
tation of satellite measurements. The asymptotic approach has been 
realized for estimating vertical ozone distribution from satellite 
measurements of backscattered solar UV radiation in the Earth’s 
atmosphere. 

L INTRODUCTION 

The backscattered solar UV radiation flux Fj in a certain 
spectral range Av centered at the wavelength n^ is formed by a 
relatively thin effective scattering atmospheric layer of the depth 
Azj. For this reason the flux Fj registered from space in the i-th 
spectral channel of a satellite instrument is connected with the 
integrated ozone content AX- in the atmospheric layer column by 
integral relation F*=f(AX i ) where f is a given integral operator. 
When measuring outgoing radiation in various spectral channels 
we obtain information about integrated ozone content in various 
atmospheric layers, i.e. about vertical ozone distribution. Such are 
heuristic grounds for remote ozonosphere sensing from space 
(Twomey, 1977). There are principal limitations for sensing accu- 
racy specified by the effective scattering layer depth: the thinner is 
the layer, the higher is the height resolution. The integral operator 
f is such that its value for given ozone profile is defined by the 
integrand function behaviour in the narrow interval Az only where 
it has evident maximum. It is well known (Olver, 1974; Nayfen, 
1981) that for such integrals are available the asymptotic represen- 
tations which connect their values with the integrand function 
values at the maximum point. In this case, other conditions being 
equal the error is specified by the interval z only. Thus, without new 
principal uncertainties being added, the remote sensing problem 
can be formulated by changing integral operators into their asymp- 
totic analogues. It allows to develop the real time algorithms for 
interpretation of the satellite measurements. This is the idea of an 
asymptotic approach to solution of the remote ozonosphere sensing 
problem (as well as any integral operator inversion problem). Let 
us attempt to realize the asymptotic idea for estimating vertical 


ozone profiles above the ozone maximum in the Earth’s atmo- 
sphere from satellite measurements of backscattered solar UV 
radiation. As a first step, suitable inverse problems are necessary to 
be formulated as a problem with a large parameter. 

2. ESTIMATING THE OZONE PROFILE AS A PROBLEM 
WITH LARGE PARAMETER 

An interaction of solar radiation with the Earth’s atmosphere 
depends on several factors (Liou, 1980; Brasseur et al., 1984). 
Apparently, the methods for inferring the vertical ozone distribu- 
tion, which should be called into action for expediting the process- 
ing of large volumes of experimental data, must be based on 
choosing of main decisive factors. From the mathematical point of 
view it may be realized by introducing the large dimensionless 
parameters with following using of asymptotic expansions. In 
particular, such method of approach in the case of horizontally 
homogeneous plane-stratified Earth’ s atmosphere (the solar zenith 
angles 0 O and the zenith angles of the receiver 0<75) allows to 
receive for values Q i (Thomas et al., 1 977) which are dimensionless 
analogues of the backscattered solar UV radiation fluxes F 
registrated from the space, the following expansion (Ishov! 
1989): 

Qi = Jo dn exp{-y i [M i X(n)/X(l) + n]}, 

where n is dimensionless height, n=n (z)=Y(z)/Y(0), n~l on 
surface at z=G, z is height, Y(z) is number of air molecules of 
atmospheric column above the height z, X(n) is total ozone content 
of an atmospheric column above the level n, X(l) is, obviously, 
total ozone content, y =ap.Y(0), a=sec0 o +sec0 is geometrical 
factor of the experiment, a [t are the i-th wavelength band 
averaged ozone absorption and molecular scattering coefficients 
respectively (Klenk, 1980, Ishov et al., 1990) and values M i of the 
dimensionless parameter M are given by 

= aoqXO)/ yi . (2) 

Accurate calculations with taking account of multiscattering 
and number of other UV solar radiation transformation processes 
m the Earth’s atmosphere (Liou, 1980) show that relation (1) is 
valid for spectral range of 250<v<300 nm with the relative error 1 % 
and can be used for interpretation of the SBUV measurements 
aboard the Nimbus-7 (Fleig et al., 1985). 
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It is easily noted that the dimensionless parameter M is equal to 
the ratio of the optical depth of atmospheric column due to ozone 
absorption to the optical depth due to molecular scattering and can 
serve a measure of mutual significance of two above attenuation 
processes of UV solar radiation in the Earth’ s atmosphere in i-th 
spectral interval. The values of parameter M for various model 
atmospheres (Ippolitov etal., 1 985) for spectral channels of SBUV 
are listed in Table 1 . One can see that the parameter M is about 2 for 
channel 302nm and continuously increases up to nearly 40 for 
channel 256nm. This fact gives ground to accounting the parameter 
M as large and the problem of estimating vertical ozone distribution 
as a problem with large parameter. 

3. THOMAS’S AND HOLLAND’S SCHEME 

Thomas and Holland have propounded a simple approach for 
expediting the processing of large volumes of experimental data for 
the purpose of evaluating the vertical ozone profiles above the 
ozone maximum (Thomas et al., 1977). Therefore, it seems natural 
at first to illustrate efficiency of the asymptotical method of 
approach to the problems of the remote atmosphere sensing within 
the framework of this scheme. 

Assuming that X(n)/X(l)=g(n) where g(n) is a given analytic 
function the inverse function G of which exists so that G(g(n))=n 
it is obtained that the equation (1) may be rewritten as 

Qi = fJdnexpj-yjfMigOi) + n]}. (3) 

Further, let T =T for all spectral channels with central wavelength 

v*?, i=l,2 L, where is the ozone absorption optical thickness 

of an atmospheric layer with upper level at the top of the atmo- 
sphere. The lower level of this layer, it is obvious, will vary for 
various spectral channels. Having chosen L different channels, we 
have fixed L different layers of the atmosphere. Apparently, the 
problem of estimating vertical ozone distribution will be solved if 
lower levels z i (or n i ) and values X(z^ (or X(n^) have been 
indicated. 

As following from our notation I i =y i M i g(n) and we have: 
n i =G(T/y i M i ), (4) 

X(n i ) = y i g(n i )/aa i ^T/acq. (5) 

Finally, it is necessary to determine the values of dimensionless 
parameter M only corresponding to experimentally obtained val- 
ues Qj for all L spectral channels. 

4. ASYMPTOTIC ANALOGUE OF INTEGRAL 
EQUATION (3) AND ITS INVERSION 

Recently, exact analytical methods of inversion of integral 
equation (3) are absent. So that we will use its asymptotic analogue 
for M— **>. From the mathematical point of view, this is the idea of 
the asymptotic method of approach to inverse problem. The 
asymptotic analogue is significantly simpler than initial equation 
since it includes linear combination of elementary functions with 
respect to M only. The error due to the substitution will approximate 


to zero when M— There are also foundations to suppose that, 
when using finite values of M in practical problem, the error will be 
not more than initial incertainty originally inherent in the inverse 
problem. The thing is half-width of integrand of (3) (resolution with 
respect to height in given inverse problem) as well as error of 
asymptotic series expansion are determinated by value of M only. 

Making substitution of integrating variable in equation (3) 
t=y i g(n), finally, we have: 

y .80) 

Qj = JdtHWexiK-M.t), (6) 

O 

where H(t) = G(t/y J )exp(-t/y i )G l (t/ y { )f y,. 


Table 1 . The values M i (2) for large dimensionless parameter M: 
1 -polar region, 2-middle latitudes winter, 3-middle 
latitudes summer, 4-tropics. 




1 

2 

3 

4 

1 

255.7 

41.45 

44.55 

40.04 

31.63 

2 

273.6 

30.23 

32.47 

29.21 

23.08 

3 

283.1 

16.51 

17.65 

16.46 

12.56 

4 

287.7 

10.58 

11.24 

10.52 

8.01 

5 

292.3 

6.47 

6.89 

6.45 

4.91 

6 

297.6 

3.46 

3.67 

3.44 

2.62 

7 

302.0 

2.05 

2.17 

2.04 

1.56 

8 

312.6 

0.52 

0.54 

0.51 

0.39 

9 

317.6 

0.34 

0.36 

0.34 

0.26 


Table 2. Estimates ofz i (km) and X^z^cm 2 ) for middle latitudes 
summer by the asymptotical method: 1 -asymptotical, 
2-exact, 3-error (%). 


No. 

z i 

1 

2 

3 

1 

46.1 

.681 (+17) 

.669(+17) 

mm 

2 

43.3 

.126(+18) 

.125(+18) 

-EK 

3 

39.9 

.270(+18) 

.272(+18) 

-1.0 

4 

37.6 

.453(+18) 


2.2 

5 

35.0 

.793(+18) 

.831(+18) 

4.6 

6 

30.8 

.161(+19) 

,172(+19) 

6.6 

7 

25.2 

.290(+19) 



8 

8.7 

. 1 35(+20) 


-64.0 

9 

5.6 

.216(+20) 

.845(+19) 


z m =24 km, 6=0.577 




Limiting oneself to consideration such functions g(n) only for 
which at t — >0 next expansion takes place 

H(t)~£c k t q ‘,-l<q 0 <q ] ..., (7) 

ko 

and using Laplas’ s method we can obtain following asymptotic at 
M— expansion instead of integral relation (3): 
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Q~Z c k r (qk + i)M qk 1 (8) 

k=o 

where T(x) is Gamma function. On the basis of (8) we can define 
M as function of Q by asymptotic at M— iteration method (Olver, 
1974) using the next iteration function 

M (n+I) ~M (0) [l + S n (M (n) )] 1/<q ° +1) (9) 

with the initial approximation for M given by 

M (o) =[cor(q 0 + l)/Q] 1/<q ° +1) (10) 

where 

s„00= ic k r(q k + 1)x< q °-^/c 0 r(q 0 + 1 ). 

k = l 


level n (Green, 1964). It can be easily shown, in this case, that 

G(t) = (st/(l - t)) 8 ,G' (t) = 6s s t 5 ‘(1 - 1)-*" 1 , 

H(t) ~ 5s 8 t 8 1 {l - (si) 8 / yf" 1 + 1(8 + l)y i +...}/ yf 

ast — » 0, so that (13) 

c 0 =8(s/ yi ) 8 ,c, = 8s 28 /y? 81 > c 3 = 6(8 + l)s 8 / 

/yf +1 q« = 5 - i.qi = 28 - i,q 3 = 5. 

By use equation (11) we obtain, finally, the next calculation 
formula for determining n^ 

n, = [T 8 'Qi /r(8 + 1)J{1 + yi r(28 + 1)/ 2r 2 (5 + 1) + 

+6(T-8-l)/s[Qj / r(6 + l)] 1/8 }. (14) 


Then stopping on the second iteration and limiting oneself to the 
three first terms in the expressions obtained, we have followed 
calculation formula: 

Mi = M \°Hl + c,r(q, + 1)[m[°> ]q 0 - q, / c 0 r(q o + 1) + 
+c 2 r(q 2 + «[M}°> jq 0 - q 2 /c o r(q 0 + 1)) (1 !) 


At last, it is necessary to estimate values of the parameters 8 and 
s in (12) and to give T, having the experimentally obtained values 
Q i only. Within shortwavelength range of spectrum, when param- 
eter M is most large, we limit oneself to the principal term of 
asymptotic expansion (8) only. In consequence we obtain the next 
asymptotic relation: 


which will be utilized below for expressing Mj through Q|. Thus, 
the problem of determination of values M i for large dimensionless 
parameter M, corresponding experimentally obtained values Q| in 
the framework of the asymptotic method, is solved. 

Table 3. The errors of the asymptotical estimates of Mj and Zj for 
tropics. 1 -asymptotical M i (exact in Table 1), 2-error 

(%), 3-exact z Jf 4-asymptotical z>, 5-error (%). 


No. 

1 

2 

3 

4 

5 

1 

31.63 

0.00 

47.17 

47.29 

-0.26 

2 

23.08 

0.02 

44.32 

44.38 

-0.16 

3 

12.55 

0.08 

40.90 

40.81 

-0.22 

4 

7.99 

0.26 

38.57 

38.30 

0.70 

5 

4.86 

1.00 

35.71 

35.37 

0.95 

6 

2.36 

9.90 

31.83 

30.76 

3.40 

7 

1.17 

25.0 

28.44 

23.93 

16.0 


5. ASYMPTOTIC ESTIMATING THE VERTICAL 
OZONE DISTRIBUTION 

Let us take as g(n) function, given by 

g(n) = n 1,8 /(n 1/8 + s), (12) 

which accords with a complete Green’s representation of inte- 
grated ozone content of vertical atmospheric column above the 


InQ ~ In T(8 + 1) / X 8 (1) - 5 In aa (15) 


which is satisfied more accurately for shortwavelength channels. 
Following (Thomas et al., 1 979) let us determine 8 from inclination 
of approximately linear dependence InQ of ln(aa) in 
shortwavelength range. The estimations show that the values 5 are 
contained between 0.5 and 0.7, so that 8 + 1 is near by position of 
Gamma function minimum. So that the values ^ obtained by 
formula (14) will be not sensitive to the errors of value 8. Parameter 
s is excluded from calculation formulas if let 

T = 8+ 1. (16) 

The results of estimating the vertical ozone distribution ob- 
tained with the asymptotic method of approach for tropics and 
middle latitudes summer are presented in Table 2 and Figure. 
Apparently, the method allows estimating X(z^ with relative error 
not more than 20% in the upper ozonosphere stretched down to the 
ozone concentration maximum height z . It appears that the values 
X(Zj) obtained with the asymptotic method near ozone concentra- 
tion maximum height are systematically underestimated, underes- 
timation incessandy increasing as approximating from top to z m . It 
is due to asymptotic character of procedure (11) used for 
determinating values of dimensionless parameter M . Analyzing 

Table 3, one can receive this result. In Table 3 are fisted relative 
errors of values Mj calculated by (1 1) and values Zj calculated by 
( 1 4) in the most unfavourable case of tropic atmosphere . It denotes 
that the relative error in the case of sought for ozone profile depend 
on value M i only. Therefore it can be easily correct, that allows 
estimating the vertical ozone distribution some below of the ozone 
concentration maximum. 
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Fig. The integrated ozone content is as the asymptotic solution of 
integral equation (3) or (1): A - g(n) = n , o - g(n)=n 1/ ®/ 
(n^+l), solid curve - model. 

5. CONCLUDING REMARKS 

Asymptotic mathematical technique proposed above permits 
using as g(n) other distinguished from equation (12) functions 
which are able more adequately describe a state of ozonosphere. In 
the particular case of g(n)=n ]/ ° we have the result which specifies 
calculation formula for obtained by Thomas et al. Although in 
this case the asymptotic method of approach does not give, in the 
main, new results, it give correct mathematical interpretation of 
them and increases the calculation accuracy. 

Other problems of remote sensing of the atmosphere also have 
an asymptotic character, so far as all of which have the same 
mathematical and physical properties. A common mathematical 
property is that integrand of the operator, which forms outgoing 
radiation, has pronounced maximum, the position of which changes 
with changing spectral channel (when scanning by wavelength) or 
direction of vision (when scanning by zenith angle and when 
tangent sounding). A common physical property is that, when 
forming outgoing radiation, registrated by satellite instrument, two 
processes are concurrent. The first is contribution to outgoing 
radiation by emitting (scattering) layer of the atmosphere and the 
second is attenuation of the outgoing radiation by upper layer on the 
path to registraring satellite instrument. A ratio of characteristic of 
these concurrent processes will play a role of large parameter in the 
problem. 
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